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Abstract 



We consider the Earth effect in the MSW analysis of the Homestake, 
Kamiokande, GALLEX, and SAGE solar neutrino experiments. Using the 
time-averaged data and assuming two-flavor oscillations, the large-angle re- 
gion of the combined fit extends to much smaller angles (to sin 2 28 ~ 0.1) 
than when the Earth effect is ignored. However, the additional constraint 
from the Kamiokande II day-night data excludes most of the parameter space 
sensitive to the Earth effect independent of astrophysical uncertainties, and 
leaves only a small large-angle region close to maximal mixing at 90% C.L. 
The nonadiabatic solution remains unaffected by the Earth effect and is still 
preferred. Both theoretical and experimental uncertainties are included in the 
analysis. 
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The recent results from SAGE, GALLEX, and Kamiokande III confirmed the deficit of 
the neutrino flux from the Sun, and thus all existing solar neutrino experiments report a 
discrepancy at some level between the neutrino flux observed and predicted If the 

experiments are correct, the implications are significant: first the standard solar model is 
excluded by the data ||[7]]. Secondly, the smaller chlorine rate relative to the Kamiokande 
rate makes astrophysical explanations difficult and excludes a wide class of nonstandard 
solar models described by a lower core temperature fU^,f|. On the other hand, among many 
proposed particle physics solutions, the Mikheyev-Smirnov-Wolfenstein (MSW) mechanism 
|T0| remains as a viable description of the data, giving a strong hint for neutrino mass. 

To explore this possibility, it is important to determine the MSW parameter space from 
the available data. The MSW effect is an experimentally verifiable hypothesis; once the 
parameter space is constrained, one can provide robust predictions for the next-generation 
detectors such as SNO, Super-Kamiokande, BOREXINO, and ICARUS. Theoretically, the 
neutrino mass and mixing angle, along with the seesaw mechanism, can be a probe of the 
physics at higher energy scales (e.g., the grand unification scale) inaccessible to laboratory 
experiments. 

If we assume two-flavor MSW oscillations into or u T , the combined data of all exper- 
iments allow two solutions [] with the squared mass difference Am 2 ~ 10~ 5 eV 2 : one in the 
nonadiabatic region with mixing angle sin 2 29 ~ 0.01 and the other in the large-angle region 
with sin 2 29 ~ 0.8, the former solution giving the better fit || (see also |5l,|lT|dT4|j). There 



is one complication in the large-angle MSW solution. For the 8 B neutrino energy range, 
the MSW parameters satisfy the resonance condition for the electron density corresponding 
to the core and mantle of the Earth. The z^'s, into which electron neutrinos convert in 



1 There is also a nonadiabatic solution for oscillations into sterile neutrinos, which is not affected 
significantly by the Earth effect. There is no large-angle solution for sterile neutrinos at 90% C.L. 
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the Sun, can oscillate back to electron neutrinos when going through the Earth. In such 
a case, the neutrino signals should be enhanced during the nighttime and can change the 
allowed parameter space of the combined fit. In this paper, by including the Earth effect, we 
improve the global analysis of Bludman, Hata, Kennedy, and Langacker, who incorporated 
theoretical and experimental uncertainties in a x 2 MSW analysis ||. 

Neutrino oscillations in the Earth have been discussed in detail by a number of authors 
1 15 1 - The survival probability P s of an electron neutrino reaching the Earth from the core 
of the Sun where it was created is calculated as in Ref. [§]; the phase difference of the two 
neutrino mass eigenstates at the entrance of the Earth can be averaged when finite energy 



bins and source broadening are included [J5J . Then the v e survival probability detected after 
propagating through the Earth is 



P s \a\ 2 + (1 - P s )\b\ 2 - {P s - -)t8Ln26{ab* + a*b) 



where a and b are the components of the unitary matrix which describes the Earth effect: 



V x (t) 



detector 
2 „;„2 



V 



a b 
-b* a* 



surface 

They depend on the energy, Am 2 , sin 2 29, the local electron density, and the path length. 
In Refs. |]15|| , the transition matrix is calculated by solving the propagation equation numer- 
ically. Here, however, we obtain an analytic formula by approximating the density profile 
with a series of either five step functions or five linear functions with respect to the Earth 



radius ||16|| . This sudden approximation is justified since the dominant effect of neutrino os- 
cillations in the Earth is the matter oscillation length becoming comparable to the neutrino 
path length, rather than a level crossing between the two energy states as in the oscillations 
in the Sun. Our calculations are in agreement with those of Baltz and Weneser and of 



Carlson |T5 | . 

The survival probability depends sensitively on the neutrino path length and thus varies 
with time, date, and detector latitude. If the parameters are in the range Am 2 /E = 1CT 7 — 
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10" 6 eV 2 /MeV and sin 29 > 10 2 , the solar neutrino signals should show not only different 
counting rates between day and night, but also a distortion of the energy spectrum with time, 
and seasonal variations during nighttime, which provides an opportunity for direct-counting 
detectors to constrain the parameter space precisely To compare with the time-averaged 
experimental results, we integrate the survival probability over each night and over a year, 
taking about thirty points for each. 

We use two data sets to constrain Am 2 and sin 2 29. One is the time-averaged rates of the 
Homestake, the combined Kamiokande II and III, and the gallium results from SAGE and 
GALLEX , which are listed in Table []} The additional information from the Kamiokande 
day-night data further constrains the parameters. We use the Bahcall-Pinsonneault solar 
model including the helium diffusion effect ||17|| . 

In fitting the data by the x 2 procedure, the experimental errors as well as all theoretical 
uncertainties due to the initial flux and the detector cross sections are incorporated. The 
X 2 value for the day-night averaged data is calculated for each Am 2 and sin 2 29 by 

X 2 (Am 2 , sin 2 29) = ]T (V-%(R° xp - R™ sw ){Rf p - Rf sw ), 

i,j=Kam,Cl,Ga 

where R^ xp and Rf ISW are the rates of the experiments and the MSW predictions for detector 
i (= Kamiokande, CI, and Ga), and V is the 3x3 error matrix. The diagonal elements 
are quadratic sums of the uncertainties: the experimental errors, detector cross section 
uncertainties (A<7j/crj =3.3% and 4% for the chlorine and gallium detectors, respectively 
fj~8||), and the flux uncertainties according to the Bahcall-Pinsonneault model ||17|| ; 



v u = ( Ai?r p ) 2 + f ^r? isw ) 2 + yi (^ R ^ Msw \ 2 

' n=pp, 7 Be, 8 B,- V <Pn J 



where R™' MSW is the MSW predictions for the n-th flux <p n (n = pp, 7 Be, 8 B, pep, hep and 



2 In principle the two predicted gallium rates differ because of their detector latitudes. However, 
the experimental uncertainties are still large and the difference is insignificant. We combine the 
two results and use 42.4°N as the latitude. 
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CNO). The theoretical flux uncertainty A0 n includes contributions from the uncertainty in 
the core temperature and in the production cross sections, added quadratically. The off- 
diagonal elements describe the correlations of the flux uncertainties among the experiments; 

n=ppjBe, s B,- \ ™ n / 

We emphasize the importance of those correlations since the dominant theoretical uncer- 
tainty from the 8 B flux (14% at la) is strongly correlated between the Kamiokande rate 
and the chlorine rate, and also because the uncertainties in the different flux components 
are correlated. Ignoring the correlation, we obtained significantly larger allowed regions, 
especially for the large angle solution 0. 

The confidence level (C.L.) contours in the Am 2 and sin 2 29 space are defined by 

X 2 (Am\sm 2 26) = xLn + &X 2 , 

where Xmin * s the global minimum of the \ 2 value in the entire Am 2 — sin 2 28 plane and 
A% 2 = 4.6, 6.0, and 9.2 for 90, 95, and 99% C.L., respectively. Caution is in order in 
determining the confidence level contours. The \ 2 procedure assumes Gaussian distributions 
for the two parameters, which is only an approximation in the present case. Furthermore the 
statistical interpretation is not clear when there is more than one allowed region. Therefore 
the confidence level contours should be considered as qualitative. 

The combined fit of the time-averaged Homestake, Kamiokande II and III, SAGE, and 
GALLEX results are displayed in Fig. [I]. For the large-angle solution the u e regeneration 
in the Earth for the 8 B neutrinos enhances the nighttime rates of the Kamiokande and 
Homestake experiments for Am 2 = 10 -6 — 10 _5 eV 2 , and enlarges the allowed region of the 
combined fit down to sin 2 29 ~ 0.1 at 90% C.L. In particular the Cabibbo-like solution 
sin 2 26 e/1 = 0.18 is now possible for Am 2 ~ 4 x 10~ 6 eV 2 . 
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The authors of Ref. [14] show larger allowed MSW regions compared to ours, because of the 



omission of the correlations of the flux uncertainties. Moreover, the uncertainties used in Ref. [14] 
are larger than the estimations of Bahcall-Pinsonneault. 



An additional constraint comes from the Kamiokande II experiment, which presented 
the results with six bins according to the angle 5 sun between the vector pointing from the 
detector to the Sun and the nadir at the detector flI5| . The first bin (cos5 sun < 0) is the 
daytime rate; the nighttime rates are for cos<5 snn = - 0.2, 0.2 - 0.4, 0.4 - 0.6, 0.6 - 0.8, 
and 0.8 - 1.0, which correspond to binning the neutrinos with different path lengths in the 
Earth. 

The Kamiokande II day-night data |19|] are consistent with no enhancement of the signal 



during the night and exclude the parameter space shown in Fig. |] [pLgf] , which is in good 
agreement with the excluded region presented in ||19|| . The exclusion comes from the com- 



parison of the six cosS sun bins, and is independent of the absolute 8 B flux. It is, however, 
subject to the uncertainties from the Earth density profile in the MSW calculations. 

In the combined fit including the Kamiokande day-night data, we treat the six cosS sun 
data as independent data points, along with the Kamiokande III, Homestake, and the gallium 
results. The normalized day-night data were taken from Ref. |19[]. They were scaled to give 
the average Kamiokande II rate; also, we included the overall systematic uncertainty (15%) 
from the energy calibration, the angular resolution, and the event selection. We checked for 
consistency by combining the six bins and reproducing the time-averaged Kamiokande II 
data (0.47±0.08 relative to the Bahcall-Pinsonneault prediction). In the x 2 fits, we calculate 
a 9x9 error matrix for each Am 2 and sin 2 29. The Kamiokande II systematic uncertainty 
is correlated among the angular bins []; the flux and cross section uncertainties are also 
included and correlated between the nine constraints. 

The result of the combined fit including the Kamiokande II day-night data is shown in 
Fig. |2|. The absence of enhancement in the nighttime Kamiokande data constrains the large- 
angle solution, leaving only a small allowed region at 90% C.L. The nonadiabatic solution 
is unchanged by the Earth effect and is still preferred. The minimum \ 2 values are 3.1 and 



We ignore possible correlations of the systematic uncertainty in the Kamiokande II and III data. 
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6.8 for the nonadiabatic and large-angle regions, respectively; both are nominally good fits 
for 7 degrees of freedom (= 9 data - 2 parameters), although the large-angle solution gives 
a poor fit without the Kamiokande II day-night data The very small x 2 contribution 
(3.0) from the six Kamiokande II angular bins is responsible for the improvement. The best 
fit parameters with x 2 value are summarized in Table [TT]. The allowed regions for 90, 95, 
and 99% CL. are shown in Fig. [3]. 

The uncertainties in the Earth density profile are not included in the above calculations. 
To examine the effect, we varied the core density and the core-mantle boundary by ±10%. 
No difference was observed for the combined fit in Fig. 
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TABLES 

TABLE I. The time-averaged solar neutrino rate relative to the standard solar model predic- 
tions by Bahcall and Pinsonneault [T^| . The Homestake result Jl| is updated from 0.26±0.04 used 
in Ref. ||. Also listed are the latitudes of the solar neutrino detectors. 





Homestake Kam II Kam III 


SAGE GALLEX 


Rate (SSM) 


0.28±0.03 0.47±0.08 0.56±0.09 

0.50±0.07 


0.44ig;^ 0.63±0.16 
0.54±0.11 


Latitude 


44.3° N 36.4° N 


43.5° N 42.4° N 


TABLE II. The best fit parameters and x 2 values for the two regions allowed by the combined 
MSW fit including the Earth effect. The Kamiokande II day-night data of the six angular bins are 
incorporated in the fit. The time-averaged data are used for Homestake, Kamiokande III, SAGE, 
and GALLEX. 




Non-adiabatic 


Large-angle 


sin 2 29 
Am 2 (eV 2 ) 
X 2 (7 d.o.f.) 


7.8 x 1(T 3 
0.55 x 1(T 5 
3.1 


0.68 
1.0 x 10~ 5 
6.8 
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FIGURES 

FIG. 1. The allowed regions including the Earth effect of the Kamiokande, Homestake, and 
gallium experiments, and the combined fit at 90% C.L. The time-averaged data are used. The 
Bahcall-Pinsonneault model and uncertainties are assumed in the calculations. Also shown is the 
large-angle solution when the Earth effect is ignored. 

FIG. 2. The allowed regions including the Earth effect of the Kamiokande, Homestake, and 
gallium experiments, and the combined fit at 90% C.L. Kamiokande II day-night data with six 
data points are included. Also shown is the excluded region (grey) using only the Kamiokande II 
day-night data. This exclusion comes from the comparison of the six angular bins and is indepen- 
dent of the 8 B flux uncertainty. 

FIG. 3. The allowed regions of the combined fit including the Kamiokande II day-night data 
at 90%, 95%, and 99% C.L. 
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